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-H . ABSTRACT 

^ ■ We cross-matched the six-year Swift/Bmst Alert Telescope (BAT) survey of Active Galactic 

I I Nuclei (AGN) with the AT20G radio survey of the southern sky, which is one of the largest 

O , high-frequency radio surveys available. With these data we investigated the possible corre- 

lation between the radio and the X-ray emission at the highest radio and X-ray frequencies. 
^ I We found 37 AGN with a high probability of association (> 80%), among which 19 are local 

Seyfert galaxies (with median redshift z = 0.03) and 18 blazars. We found that w 20% of the 
AGN detected in hard X-rays are also bright radio sources at 20 GHz, but the apparent correla- 
tion between the radio and hard X-ray luminosity is completely driven by the different median 
K*" ■ redshifts of the two subgroups of AGN. When we consider only the local Seyfert sample we 

OO ' find no evidence of a correlation between their 20 GHz and 15-55 keV power Therefore it 

appears that at high frequencies the radio-X connection, which had been previously observed 
at lower frequencies, disappears. The disappearance of the radio-X correlation at high radio 
and X-ray frequencies could be tested through Very Long Baseline Interferometry and the 
■^ . use of the Nuclear Spectroscopic Telescope Array (NuSTAR) satellite. 

ff^ ■ Key words: radiation mechanisms: non-thermal, quasars: general, radio continuum: general. 

X-rays: general 

>: 

X 

3 1 INTRODUCTION The most accepted scenario for compto nization requires a popula- 

" " " tion of hot electrons (i.e. a corona, see iHaardt & Maraschilll99ll) 

Active Galactic Nuclei (AGN) are among the most powerful ,i . ii j ■ i .. ,i ■ ^ ,^ r -^ ^- i j-- r> 

" '^ that sandwiches the inner part (tens oi gravitational radii it„ — 

sources of radiation across the whole electromagnetic spectrum, ^>^;2x^i ■ ii ii-i-i , i 

GM/c ) of the accretion disk, and which is kept hot by magnetic 

from 7-rays (up to TeV energies in some cases) down to radio . »,,-,, r, ■ , , ■, 

reconnection. Additionally, an outflow or a jet could contribute 

waves. Accretion onto a central supermassive black hole (SMBH) , ,i_ » » i -.» j ■ .i . j- r 

^ to the total emitted power, more prominently at radio frequencies 

IS generally accepted to be one of the dominant physical mecha- , , , ... , , , ,, „ 

. through sync hrotron radiation and less markedly at 7-ray irequen- 

msms that power these systems. Another source of energy is be- , rr riu„, ,1 r; : — rTT^TTtTl 

,. , , , . , ^ , ^, ,„,, . ,^ , . , . cies (see e.g. lTeng et al.ll201 IklLenain et alJ201(J) . 
lieved to be the rotational energy of the SMBH itself, which in 

some cases can exceed the gravitational energy extracted through when AGN also produce relativistic jets which are pointing in 
accretion. tjje direction of the observer, they are referred to as blazars. In this 
The radia tion output of quasars and Seyfert galaxies peaks in subclass of AGN the SED shows two clear peaks, or "humps". In 
the ultraviolet (jShakura & Sunyaev||l973|), producing the so-called leptonic models, all the radiation we see coming from the central 
"big blue bump" in the spectral energy distribution (SED), while object is generally thought to be generated by a single population of 
only a fraction of the energy available through accretion is repro- electrons, i.e. directly or indirectly by the jet itself. This is achieved 
cessed into high energy emission (from keV energies to hundreds via synchrotron emission producing the lower-energy hump and 
of keV) via thermal comptonization, which produces a second, via a combination of external compton (EC) and Synchrotron-self- 
less dominant peak in the SED, and which is generally thought t o Compton (SSC) producing the high energy hump. Only in a hand- 
cause the cosmic X-ray background (CXB, [Giacconi et al.||T96a). fui of cases the thermal radiation coming from the accretion disk 

is powerful enough to "pierce" through the synchrotron hump, and 

could be used to estimate the mass and accretion rate of the AGN. 

* E-mail: dburlon@physics.usyd.edu.au All blazars are powerful radio sources, and the relative contribu- 
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tion of the two peaks in the SED as well as the position of their 
synchr otron peaks classify them into subclasses of radio (and total) 
power ( lFossatietal.lll998h . 

In practice, a popular method of classifying AGN uses the ra- 
tio of the radio (at 5 GHz) to optical (at 4400 A) fluxes, through 
the so-called radio-loudness parameter R = Fradio/FopUcai- Ini- 
tially, the distinction between radio-loud (R > 10) and radio- 
quiet (J?<10) objects appeared sharp (Stri ttmatteret al.l Il980l : 
iKellermann et alj|l989l . Il994) if derived from the distribution of 
the R parameter of AGN present in the Palomar Bright Quasar Sur- 
vey. Up to five other similar prox ies of R have bee n defined in the 
literature (for a recent review see Elvis et al.ll2012l and references 
therein). The existence of a bimod ality has been ever s ince debated 
in the literature in favor ( s ee e.g. Miller et al.lll990t Stocke et al.l 
1 19921 : iHooper et al.l 1 19951: llvezic et alJ |2002| : Jjiang et alJ l2007l) 



interpretation was discuss ed in relation to the scale invariance of 
black hole jets (see also H einz & Sunvaevl2003h. Other authors find 
a slope close to unity (Canosa et al. 1999: Brinkmann et al. 2000; 
Salvato et al. 2004; Panessa et al. 2007) or a completely different 



slope Lx oc L?j® jWang et alj 200q) similar to what observed in 



or against it (see e.g. Condon et al. 1 980l : I White et al.l I200C ; 
ICirasuolo et al . 2003; Rafter et al. 2009). A possible source of con- 
fusion is the different size, sensitivity limits, cuts in redshift, optical 
and radio frequencies used to classify AGN through the R parame- 
ter. 

Recently iMahonv et alj ( 120121 ) used a soft X-ray selected 
sample of several hundreds quasars and did not find any evidence of 
bimodality. Seyfert galax ies are traditionall y considered radio-quiet 
sources but, for example. IHo & Pend ( l2001i) showed that consider- 
ing their nuclear luminosities, most Seyfert galaxies are radio-loud. 

To account for this dichotomy, it was proposed that just rapidly 
spinning black holes can produce relativist ic jets responsible for 
non-thermal emission (e.g. iBlandfordll 199(1) . Nonetheless, in later 
years the evidence for high spin came from a radio-qu iet object 
(MCG-6-30-15, Tanaka et al. 1995). Other works (Na ndra et al.l 
120071 : [Patrick et al.ll2012i) showed that broad relativistic lines are 
frequently observed in Seyfert galaxies and in a handful of cases 
(e.g. NGC 3516, NGC 7469) they could also constrain the spin. 
Other results are at odds with the idea of radio-quiet objects host- 
ing slowly spinning black holes. One example is the high ef- 
ficiency required b y accretion in order to account for the inte- 
grated CXB CElvis et al.l 120021) . which is more easily explained 
if the BH are spinning. Other examples are X-ray observations 
that indicate that the emission extends well within the innermost 
stable orbit of non- rotating black holes (e.g. Fabian etal. 2002; 
llwasawa et al.lll996r) . Moreover, both Seyfert galaxies and low- 
luminosity AGN (LLAG N) were shown to be radio emitters at 
some level. For instance, IHo & Ulvestag ( 1200 ih have shown that 
85% of Seyfert AGN are detected at 5 GHz, with a wide range 
of radio powers and morphologies, where structures could be re- 
solved by very long bas eline interferometry (VLBI). J Nagar et al.l 
( |2002|) . lUlvestaj ( l2003h and lOiroletti & Panessal ( |2009|) showed 
clear structures (blobs or knots) at the sub-parsec level. The origin 
of radio emission at this level is not yet fully understood. In most 
cases it is associated to synchrotron emission just as in the r adio- 
loud case dGiroletti et al.l2005l : llshibashi & CourvoisienBOlOl) and, 
if so, should highlight the presence of a jet or outflow also in this 
class of objects. Alternatively, thermal free-free emission is viable 
option (see e.g. iGallimore et al.ll2004l for the radio emission of 
NGC 1068, among others). 

The radio-X connection has been addressed by different au- 
thors, but given the different energy ranges used in the high-energy 
band coupled with different radio frequencies, the results are of- 
ten discordant (and will be discussed in S ect. 4). For instance, in 
the seminal work by iMerloni et al.l ( I2003B a relation between the 
radio luminosity at 5 GHz and the X-ray 2-10 keV luminosity 
Lr oc Lx^ (i.e. Lx oc L]f) was found and a solid theoretical 



few Galact ic black holes I Corbel et al.ll2013h . While completing 
this paper, ( IBonchi et alj|2013[) reported an updated version of the 



Fundamental Plane, with a slope even harder than the original one: 

r ^ rO.39 
Lr oc Lx ■ 

The aim of this work is to study for the first time the radio- 
X connection in local Seyfert galaxies in the hard-X ray, high- 
frequency radio regime. We want to test whether the correlation 
between the radio and X-ray power at high radio and X-ray fre- 
quencies is as strong as at lower frequencies, or conversely tends to 
vanish (as it seems to be the case when high resolution Very Long 
Baseline Interferometry observations are used - Panessa & Giro- 
letti, 2013, MNRAS submitted). To these aims we cross-match the 
sources present in the BAT survey (15-55 keV) with the sources de- 
tected by the high frequency (20 GHz) Australia Telescope Com- 
pact Arra y (ATCA) survey. A similar complementary work is in 
progress jPanessa et al.ll201 11) adopting the sample of AGN de- 
tected by the INTEGRAL satellite, and the 1.4 GHz radio lumi- 
nosity. Although our primary focus is on Seyferts, we also con- 
sider the blazars found through the cross-matching, in other to in- 
vestigate the possible gap in the radio-to-gamm a ray connection 
in blazars recently studied JMahonv et al.l uOlOl : iGhirlanda et al.l 
l20IOl : lAckermann et al.ll201lll 

Another approach to the present work is to investigate whether 
the correlations arise when long-term averages are used, or con- 
versely if there is a radio-X connection at all scales. When at lower 
radio frequencies the "core" luminosities (albeit the "core" could 
mean jet structures of the order of kpc) are confronted to the high- 
energies the correlation holds both for pointed and survey (i.e. av- 
eraged) X-ray observations. Choosing higher radio-frequencies al- 
low us to be less contaminated by the extended structures and there- 
fore a shorter (but still somewhat averaged) radio output. We ex- 
pect that if we compared hard X-rays fro m a survey telescope with 
lower radio frequencies (see an example in lPanessa et al.l201 ih . we 
would still recover a strong correlation, possibly because we would 
be confronting longer term averages, instrumental/observational 
in the case of X-rays and physical in the case of radio waves. 
Nevertheless, we note that there is a still large difference in the 
time scale for X-ray emission (years, for all-sky surveys) and ra- 
dio waves (several thousand years in kpc-scale regions). The best 
test we envisage for the future will be to investigate the radio-X 
connection using foc ussing hard X-ray telescopes like NuSTAR 
I IHarrison et al.ll2010l) and VLBI radio luminosities, which should 
both probe at the same time the very inner part of the AGN. 

In Sect. 2 we briefly discuss the two surveys and the cross- 
correlation technique. We present our results in Sect. 3 and discuss 
them in Sect. 4. We summarise our findings in Sect. 5. A standard 
flat Universe with h = SIa =0.7 is assumed. 



2 THE SAMPLE 

We give a brief description of the two blind surveys we used for 
this work, but we refer the reader to the respective original pub- 
lications for a detailed discussion of both surveys' strategies and 
peculiarities. In §[23] we describe the cross-correlation technique 
we adopted. 
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• Matches (>95%) 

♦ ESOS49-G049 
■ ESO1S7-G023 



^'ObX,y °'' 



Figure 1. Association (posteiior) probability between the Swift/BAI and 
AT20G samples versus radio flux density and hard X-ray flux of the 
matched sources. Two 'outliers' are clearly visible in the low-flux low- 
probability section of the data cube. These are the only two AGN for which 
the association probability is between 80% and 95%. 



2. 1 The SwiftfBAT survey 

The Burst Alert Telescope (B AT; Barthelmy et al.l2005h instrument 
on-board the Swift satellite l IGehrels et al.ll2004l) is a coded mask 
instrument originally thought to alert the X-ray and Optical-UV 
cameras to follow-up 7-ray burst afterglows. While surveying the 
sky waiting for a trigger, the BAT continuously monitors up to 
80% of the sky thanks to its wide field of view (120° x 90° par- 
tially coded). After several years of deep exposure, the BAT sur- 
vey covers almost uniformly the whole sky down to the limiting 
sensitivity. In our analysis we used the data accrued in 6 years 
of exposure, down to a sensitivity limit ~ 6 x 10~^ ^ erg/cm^/s. 
The BAT survey is presently (see lAiello et al.ll2012l for the lat- 
est release) the deepest blind scan in the sky above 10 keV mak- 
ing it the optimal candidate for the goal of confronting the hard 
X-ray and radio properties of AG N. Data screening a nd process- 
ing were describe d in details in ( lAiello et al.l I2OO8I A08a) and 
dAiello et alj|2008l . A08b) but we summarise hereafter the main 
characteristics of the survey. The chosen energy interval for (con- 
servatively) computing the fluxes was 15-55 keV, while the spec- 
tra are traditionally extracted in the 15-200 keV energy range. 
The all-sky image was obtained by averaging the (weighted) per- 
pointings, resulting in a final sample of ~ 350 AGN, comprising 
mostly Seyfert-like galaxies (284) and blazars. Just ~ 5% of the 
sources of the BAT catalogue still escape identification, and this 
incompleteness is independent of the Galactic latitude (see e.g. 
lAiello et al.l 120121) . Each of these sources has a signal-to-noise 
ratio SNR>5a, and was identified by cross-matching it aga inst 
the previous catalogues of Tuelleretal. (2008), Cusumano e t al.l 
( |2010|) . IVoss & Aiellol ( l2010h . and iBurlon et all l l201ll). When- 
ever possible , optic al identifications were taken fr o m | Masetti et al.l 
( I2OO8I |2009|. I2OI0I) as described in lAiello et al.l ( I2OI2I) . For the 
blazar s we checked their c lassification against the BZCATlj cata- 
logue JMassaro et al]|2009l) . 
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Figure 2. Scatter plot of the radio spectral indices of the all AT20G sources 
are shown in red crosses. Seyfert galaxies (blazars) are shown with circles 
(squares). The two panels show the (arbitrarily reseated) projected distri- 
butions along the respective axis of the whole AT20G population (red his- 
togram) and the whole cross-correlated subsample (black histogram). 



2.2 The AT20G survey 

The AT20G survey dMurphv et al.ll2010l) is one of the largest blind 
surveys at high radio frequencies ever conducted from ground. It 
covers ~ 20, 000 deg^ to a limiting flux density of 40 mjy/beam. 
The key feature of the survey was its two-phase strategy, where a 
fast blind scan (possible due to the fast-scanning mode of the Aus- 
tralian Telescope Compact Array) was followed by regular snap- 
shots by the ATCA. 

In the first phase the ATCA scanned the whole southern sky at 
20 GHz, with a speed of 15° min^^ in declination at the meridian. 
The scanning strategy exploited the Earth rotation, while sweeping 
regions of the sky 10° — 15° wide in declination. The rms sensitiv- 
ity of the first phase was ~ 10 mJy. In the second phase, the sources 
brighter than 5a where followed-up with an hybrid array configu- 
ration and two 128 MHz bands, combined during data processing 
at the central frequency of 19.904 GHz (that we refer to as 20 GHz 
throughout the paper). The sources with declination < —15° were 
observed also at 4.8 GHz and 8.64 GHz (which we refer to as 
5 GHz and 8 GHz), with an east-west configuration. The exclu- 
sion of sources with declination > —15° for the lower frequencies 
is motivated by the poor {u, v) coverage of east-west arrays near 
the equator. 



2.3 Source Associations 

For the positional cross-matching of the BAT and AT20G source 
catalogs we adopt a Bayesian method already used in cross- 
match ing the Fermi AGN with the AT20G survey dGhirlanda et al.l 
UOIOJ) and first dev eloped for the identification of the Fermi sources 
dAbdo et al.ll201(l and references therein). This method, given the 
BAT source position and accuracy (which we assume to be half 
of the point spread function) assigns a probability to every radio 
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Figure 3. Distribution of the parent populations flux, compared witli the associated AGN. The distributions have been rescaled. (Left:) Sn'i/i'/BAT flux in the 
energy band 15-55 keV for the whole 6-year catalog (thin line) and the associated subsample (thick dashed line). (Right:) AT20G flux densities at 20 GHz, 
with the same notation as left panel. 



source found within the BAT error circle. This probability is com- 
puted accounting for the surface density of radio candidates around 
each BAT source. To this aim, the prior probability of the radio 
sources (i.e. the confidence that the radio catalogue contains the 
counterpart of the X -ray selected AGN ) is calculated as per ap- 
pendices B and C of Ab do et al.l j2010f) , and is assumed equal for 
all the sources. In practice, this translates into calculating the aver- 
age number of false associations between the AT20G catalogue and 
100 ran domly gener ated mock X-ray cata logues . Simi larly to what 
done in lAbdoetai] (2010) and Ghirlanda ltal] ( |2010|) . we call an 
'association' a unique match (i.e. a single radio source within the 
BAT error circle) with a posterior probability greater than 80%. We 
comment bel ow on some of the so urces that are flagged in the orig- 
inal paper by Murphy et al.l OOlOl) . While this method allows us to 
derive a figure of merit for the associations found in the two cat- 
alogs, we note that by only searching for positional matches (i.e. 
radio sources within the BAT error circle), we could derive a simi- 
lar -albeit smaller- number of associations. 

We find 37 BAT AGN with a radio counterpart in the AT20G 
survey. Not surprisingly the blazar class has a detection rate of 
~ 65%, which is much higher than the equivalent for Seyfert 
galaxies, i.e. ~ 14%. These are listed in Table [T] along with the 
posterior probability of the association. About half of the associ- 
ations (19/37) are classified as Seyfert galaxies, roughly equally 
distributed into the two subclasses of optically Seyfert "type 1" (7, 
plus 1 Syl Narrow Line) or "type 2" (11). There are 18 AGN classi- 
fied as either BL Lacs (1) or FSRQ (15) blazars, and 2 classified as 
generic blazar type. The absolute majority of the FSRQ class is no 
surprise, given that our parent sample is the Swi/C/BAT one, which 
is well known to pick preferentially this class of AGN. Note that all 
but two Seyfert 2 AGN (i.e. ESQ 549-G049 and ESQ 157-G023) 
have a probability in excess of 95%, and that both the less certain 
associations are close to the detection limit of both surveys. This is 
represented in Fig. [T] where the two sources are clearly displaced 
from the bulk of the associations. 

Properties of the cross-matched sample: One of the results of 
the AT20G survey is that a large fraction of sources are not well 
fit by a single power law even in the small bandpass of the ATCA 
telescope. In Fig.[2]we show the distribution of the spectral index of 
the several thousands sources of the AT20G for which the spectral 



index could be comp uted (see for comparison Murph v et al-bOld : 
iMassardi et al.ll201 ih . In the same plot we show the radio spectral 
indices of the BAT sources associated to the radio counterparts. The 
spectral indices whe re calculated in the standard way taking from 
iMurphv et alj ( 120 lOh the flux densities (5) at frequencies vi and 

_ loguAS{i'i)/S{v.2)) 
logu){vi/v2) 

Fig. |2] shows that the associations we found do no show a particu- 
lar spectrum, i.e. they populate all the four quadrants of the plane 
identified by the dashed lines corresponding to a = —0.5 which is 
often adopted as a s eparation between hard and soft radio sources. 
IChhetri et akl ( 1201 2|) pointed out that AT20G is a high frequency 
survey and hence it has a larger fraction of flat-spectrum com- 
pact radio sources (~ 80%) than the surveys at lower frequencies 
(~ 25%). 

We checked also whether the association method tended to 
pick the brighter end of the distribution of both populations. In 
Fig.[3]we show that albeit this is not true for the hard X-ray popu- 
lation, it seems that the radio part of our association picks up pref- 
erentially brighter sources. When we compare the 5w(/f/BAT flux 
distributions (left panel), it is apparent how this work is sampling 
AGN which are on average distributed exactly as the parent popu- 
lation. When we look at the 20 GHz flux distributions, the associ- 
ated AGN have a flat distribution. Note that we chose to compare 
the flux densities at 20 GHz because this is the most sampled fre- 
quency of the radio surv ey (the follow-up com pleteness is > 90% 
only above 100 mjy, see'Massardi et al. 20n)^ 



Notes on 'extended' source}^ Massardi et al.l ( 120081) reported 
that the radio galaxy Fornax A (that is not present in our sam- 
ple) is the only bright extragalactic source known to be missing 
from the AT20G. This incompleteness is a result of its 20 GHz 

^ In the original AT20G catalog 5 of our sources were originally flagged as 
'poor' . Their flux density measurements are anywa y unaffected due to the 
method (refen'ed to as Triple Product Method; see lMurphv et alj|2010l for 
additional details) used, which is robust against the effects of decorrelation 
(phase fluctuations of the visibilities). These flagged data have no signifi- 
cant effect for the analysis. 
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Table 1. Cross-matched sample. The sources are Hsted with increasing AT20G name. The 5w;yf/BAT flux is in units of erg/cm^ /s, while the radio flux densities 
are in mjy. The BAT fluxes have been de-absorbed with the recipe of Burlon et al. (201 1) where logAf^f was in excess of 23.5, and these sources are marked 
in boldface accordingly. References are in the text. 



Name 


Type 


AT20G Name 


Redshift 


Prob. 


Fluxi5_55 


S20 


O"S20 


S8 


0"S8 


S5 


0-S5 


^BAT 












[cgs] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 


[mJy] 




NGC 612 


Sy2 


J013357-362935 


0.0297 


0.994 


3.535e-ll 


440 












1.63 


NGC 1052 


Sy2 


J024 104-08 1520 


0.0050 


0.952 


1.359e-ll 


1199 


59 










1.47 


NGC 1068 


Sy2 


J024240-000046 


0.00378 


0.995 


9.985e-10 


474 


17 










2.23 


PKS 0312-77 


FSRQ 


J03 1 155-765151 


0.2230 


0.987 


9.237e-12 


1238 


81 


1292 


67 


924 


46 


2.08 


RES 0471 


Syl 


J034730-303521 


0.095 


0.979 


7.383e-12 


55 


2 


42 


2 


57 


3 


2.04 


IGR J03532-6829 


BLAZ 


J035257-683117 


0.0869 


0.987 


9.778e-12 


68 


3 


114 


4 


165 


5 


2.30 


ESQ 549-G049 


Sy2 


J0403 14- 180953 


0.0262 


0.874 


1.394e-ll 


54 


4 


106 


7 


146 


7 


1.75 


RXJ0405.5-1308 


FSRQ 


J040534-130813 


0.571 


0.960 


8.180e-12 


1576 


77 










2.14 


ESQ 157-G023 


Sy2 


J042203-562127 


0.0433 


0.800 


9.956e-12 


95 


5 


105 


5 


85 


4 


1.67 


LEDA 177818 


Sy2 


J044436-280914 


0.147 


0.983 


9.755e-12 


221 


10 


1431 


19 


2245 


46 


2.17 


Pictor A 


Syl 


J05 1949-454643 


0.0351 


0.985 


2.136e-ll 


6320 


110 










1.90 


ESQ 362-G021 


BLLAC 


J052257-362730 


0.0553 


0.983 


1.612e-ll 


3909 


156 


6566 


244 


9066 


342 


1.85 


PKS 0524-460 


FSRQ 


J052531-455754 


1.479 


0.988 


9.642e-12 


483 


24 


698 


35 


845 


42 


1.22 


SWIFT J0539.9-2839 


FSRQ 


J053954-283956 


3.103 


0.986 


1.009e-ll 


677 


44 


1143 


60 


1278 


64 


1.55 


PKS 0558-504 


SylNL 


J055947-502652 


0.1369 


0.983 


1.014e-ll 


78 


4 


70 


4 


104 


5 


2.04 


PKS 0637-752 


FSRQ 


J063546-751616 


0.6350 


0.988 


1.079e-ll 


3142 


206 


4358 


227 


4804 


240 


1.92 


3C206 


Syl 


J083950-121434 


0.1978 


0.995 


1.216e-ll 


323 


16 










2.04 


PKS 0921-213 


BLAZ 


J092338-213544 
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0.989 
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22 
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20 
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20 


2.24 


QSOB1127-145 
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1.187 


0.990 
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1.75 
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0.5361 


0.995 


1.233e-ll 


20024 
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1.50 


NGC 4945 


Sy2 


J130527-492804 
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0.995 


5.263e-10 
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36 


786 


39 


1566 


78 


1.47 


IGR J13 109-5552 


Syl 


J131043-555209 


0.1040 
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1.219e-ll 
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1 


321 


16 
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16 


1.40 


CenA 


Sy2 


JI32527-430104 


0.00182 
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5.687e-10 
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1.85 


Cir Galaxy 


Sy2 


J141309-652020 


0.00144 


0.988 


1.035e-09 


97 


5 
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7 
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15 


2.09 


PKS 1510-089 


FSRQ 


J15 1250-090558 


0.3599 


0.970 


2.656e-ll 


2933 


142 










1.36 


PKS 1549-79 


FSRQ 


J155658-791404 


0.1495 


0.976 


1.170e-ll 
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41 
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99 


3326 


166 


2.47 


SWIFT J1656.3-3302 


FSRQ 


J165616-330207 


2.400 


0.984 


2.610e-ll 


287 


14 


241 


12 


220 


11 


1.60 


PKS 1830-21 


FSRQ 


J183339-210341 


2.500 


0.987 


2.343e-ll 
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360 










1.50 


NGC 6860 


Syl 


J200848-610938 


0.01488 
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3.309e-ll 
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Sy2 


J205202-570407 


0.0113 
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4.034e-ll 
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6 


230 


12 
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23 


1.89 


QSO B2052-47 


FSRQ 
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1.488 
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1.015e-ll 
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57 
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199 


3204 


160 


2.25 


PKS 2126-158 


FSRQ 


J212912-153841 


3.280 


0.987 


1.268e-ll 
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54 
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82 


1528 


76 


1.62 


PKS 2149-306 


FSRQ 


J215155-302753 


2.345 


0.987 


3.526e-ll 
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1.48 


NGC 7213 


Syl.5 


J2209 16-47 1000 


0.0277 
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2.462e-ll 
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6 
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8 
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0.05639 
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2.112e-ll 
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1.99 


PKS 2227-08 


FSRQ 


J222940-083254 


1.561 


0.987 


9.945e-12 


3230 
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1.52 


PKS 2356-61 


Sy2 


J235844-605246 


0.0963 


0.986 


9.463e-12 


311 


12 


1482 


21 


2124 


41 


1.80 



flux coming from the lobes being resolved by the ATCA beam. In 
our sample 6 sources (2MASX J04450628-2820284, NGC 1068, 
Circinus, PKS 2356-61, IGR J03532-6829, ESO~362-G021) are 
flagged as 'extended', while 3 additional ones (Cen A, Pictor A, 
NGC 612) are additionally flagged as 'large and extended'. Note 
that the AT20G is known to be biased against extended sources 
due to a reduced sensitivity (50% reduction in amplitude for source 
si ze > 45 arcsec) and all these sources were targeted specifically 
bv lBurke-Spolaor et al.l ( 12009) . see in particular their Table 2. Their 
total integrated flux measurements are reported in the catalogue we 
are using in this work, so that when we compare with the compact 
sources, we are not underestimating the flux densities of the former. 
In practical terms, we are using the flux coming from the innermost 
1-2 kpc of the source (e.g. for Cen A the scale is 0.037 kpc/arcsec, 
which means that the 20 GHz flux densities refer to the innermost 
1.57 kpc, i.e. the core of the active galaxy and the initial part of the 
jets that extend for > 1 Mpc). Also Murphy et al. (2010) comment 
further on this, in the context of associating the radio emission with 



the core of NGC 1068 (and other close galaxies), and point out that 
the reported flux appears to be associated with the active nucleus. 



3 RESULTS 

We explored the possible correlation between the hard X-ray flux 
(15-55 keV) and the radio flux density at 20GHz. Fig.|4]shows the 
scatter plot for the associations. There is no correlation between 
the two quantities for the Seyfert galaxies as also confirmed by the 
Spearman's rho correlation coefficient Ts — 0.38 (corresponding to 
a chance coiTelation probability P — 0.11). This correlation is even 
weaker if we include the blazars. Therefore we find no evidence of 
a relation between the fluxes. 

Note that four associations are Com pton-thick sources : 
NGC 612 and NGC 1068 (both present in Isurlon et aDbOllh . 
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Figure 4. Scatter plot of the Sh'i/i'/BAT flux in the energy band 15-55 keV 
vs. 20 GHz flux density. The sensibilities of the two surveys are represented 
with the grey shaded area. Seyfert galaxies are plotted with blue open (red 
fill ed) circles if their flux from the six-year BAT survey (unabsorbed as 
per lBurlon et al.ll201 11) is used. Blazars are represented with black squares 
(FSRQ), crosses (BL Lac), and diamonds (unclassified). 



NGC 4945, and Circinus. We show their absorbed and unabsorbeqj 
X-ray fluxes (connected by a vertical dashed line) in Fig|4] 

We over-plotted in Fig.|4]the two limiting fluxes of both sur- 
veys as a grey shaded areas. The absence of points in the upper part 
of the plane does not seem affected by a simple selection effect, 
while the lower right one is likely affected by the BAT sensitivity. 

However, the study of the correlation between the power out- 
put at different frequencies is often done in the luminosity space 
rather than in the flux space. We computed the fc-corrected BAT 
luminosities via: 



I/BAT = 47rdL(2) 



Fl5- 



{1 + zY 



where T is the photon index of the individual source fitted in the 
whole BAT energy band. Note that the usual notation of the X-ray 
photon indices is such that F{u) oc u^~^ . The fc-corrected radio 
luminosity at the i-th central frequency (19.904 GHz, 8.64 GHz 
and 4.8 GHz) was obtained from the flux density ^(i^) via: 



LaT2<3 = 'iTxdh(z) 



Vi X S{Vi) 

(l + z)i+° 



where q is the energy spectral index {S(y) oc i^"), so that we are 
calculating both luminosities in units of erg/s. We chose to use the 
spectral index between 5 GHz and 20 GHz as a proxy of the radio 
spectrum in the available ATCA band. We used the median spectral 
index for the sources which lacked the flux density m easurement at 
5 GHz . A caveat for the more distant sources is that lChhetri et al.l 
( 1201 21) showed that the AT20G sources tend to show a turn over 
at the rest-frame frequency of about 30 GHz. This curvature effect 
does not significantly impact the fc-corrected luminosities for all 



Table 2. For each frequency block the first and second line are the Spear- 
man's correlation coefficient and corresponding null hypothesis probability. 
The third and fourth lines show the corresponding partial correlation analy- 
sis results (labelled with a (z)), as described in the appendix. 





All 




Seyferts 


Blazars 


20 GHz" 


(37) 




(19) 


(18) 


-^BAT-iAT20 
-^BAT-iAT20 (z) 


0.94 
7.99x10" 

0.38 
1.96x10 


19 
-2 


0.86 

1.76x10^'^ 

0.07 

0.78 


0.84 

1.11x10-5 

0.13 

0.59 


8 GHz" 


(22) 




(11) 


(11) 


-^BAT--^AT20 
-'jBAT--^AT20 (z) 


0.95 
6.74x10" 
0.23 
0.29 


12 


0.94 

1.12x10^5 

0.37 

0.25 


0.84 

1.33x10-^ 

0.07 

0.83 


5 GHz'' 


(22) 




(11) 


(11) 


LBAT-iAT20 
-'^BAT-iAT20 (z) 


0.95 
8.55x10- 
0.26 
0.23 


12 


0.94 
1.12x10-5 

0.42 
0.19 


0.84 

1.33x10-=* 

0.05 

0.87 



^ The numbers in parenthesis represent the number of AGN used to com- 
pute the statistics. 



but the sourceQ which are at high redshift (z>l) and harder than 

The result is shown in Fig. |5] The I/bat-I/AT20 coiTelation 
appears stronger with respect to the flux plane, as can be seen in the 
corresponding Spearman's rho correlation coefficients and chance 
correlation probabilities (see the top panel of Table |2)- A fit with 
the bisector method at 20 GHz yields 

logLx = 0.68(±0.11) logLn + 15(±5) 

for the Seyfert sample, which is consistent with the slightly 
flatter result we obtain including also the blazars logLx ~ 
0.75(±0.05) logLn. + 12(±2). 

However, the common dependence of the radio and X-ray lu- 
minosities with redshift should caution about the significance of 
the correlations. In order to remove the effect of z we performed 
a partial correlation analysis, as detailed in the appendix. We find 
that the correlation is not significant for the sample of Seyfert-like 
AGN (chance correlation probability is 78%). The apparent relation 
is therefore completely driven by their distance modulus. A caveat 
when considering the whole sample: from a purely statistical point 
of view the null hypothesis probability is ~ 2 x IQ-^, correspond- 
ing to a confidence of 2.4 a. Nonetheless, since the redshift distri- 
bution is almost bimodal, we are effectively computing the partial 
correlation among two different populations. The Seyfert one with 
a lower median flux and distance versus the blazar one, with a much 
higher median flux and redshift. The results of the analysis in all the 
three available radio bands is summarised in Table |2] 



■^ For NGC 4945 and Circinus we estimated the suppression factor us- 
ing the prescri ptions of J Bur lon et al. (2011), which are consistent with 
the findings of iBrightman & Nandral 120111) . To this aim we used a fidu- 
cial value of Nu take n from the most recent literature ( lYaqoob..2012. and 
iDella Ceca et alj2008l respectively) 



* For SWIFT J1656.3-3302, SWIFT J0539.9-2839, PKS 2126-158, 
QSO B2052^7, PKS 0524^60 the radio luminosity is strictly speaking 
a lower limit. 
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Figure 5. Swi/i'/BAT luminosity vs. AT20G radio luminosity at 20 GHz. Different sources as per the legend and Fig.|4] The best (bisector) fits are shown with 
a dashed blue, dash-dot, dash-triple-dot line for the Seyferts, all sample, and blazar respectively. The dotted lines represent three different Kx values of 0.01, 
0.1, 1 and are conversely expressed in the bottom left of the comer as relative radio-X luminosities. 



4 DISCUSSION 

There is a natural scenario in which all X-ray emitting AGN are 
also radio emitters at some level, i.e. a scenario in which all AGN 
are able to accelerate a jet, even though only a small fraction 
(~ 10%, the genuine radio-loud) are able to successfully launc h 
such a jet. In this "aborted-jet" scenario JGhisellini et al.l 120040 , 
most AGN are indeed unable to launch the jet due to its subrel- 
ativisic speed, which make the jet collapse back and account for 
both radio and high-energy emission. The difference with the more 
well-known accretion scenario is due to the heating mechanism for 
the electrons responsible for the X-rays, which in the case of an 
aborted jet, is predicted to be localised along the (failed) jet axis. It 
is worth nothing that both scenarios are not necessarily in competi- 
tion with each other, and may well be both contributing at different 
levels in the radio and high-energy bands. 

It is not surprising that lot of effort has been spent investi- 
gating the reciprocal contribution of radio emission and X-rays. In 
recent years numerous groups have tackled this issue from different 
perspectives. D ifferent authors used soft X-ray f Canosa et alj 19991 : 
iBrinkmann et al. 2000 : Salvato et al. 2004), and found a slope of 
the correlation which is close to unitjQ Some differences among 
the different works emerge as a consequence of the fitting method, 
and/or of the inclusion of some AGN which might be contaminated 



^ The slope is consistent with unity for IBrinkmann et al.l J2000l) if the 
radio-quiet sample is considered. 



by a starburst compo nent, and/or the defi nition of radio-loudness 
(as discussed before). ICanosa et al.l ( 119991) adopted monochromatic 
X-ray luminosi ties at 1 keV for a sample of ~ 40 (2 < 0.1) 
AGN. Brinkma nn et al.l ( 120000 chose 2 keV luminosities for ~ 60 
{z < 1) AGN , respectively. These were compared to the 5 GHz 
luminosity and a significant correlation was found (even if the 
ef fect of distance migh t have been underestimated). In the work 
of ISalvato et"ZI ^2004^, 53 local {z < 0.1) AGN were cross- 
correlated in the X-ray energy band of 0.5-2 keV versus 1.5 GHz 
monochromatic luminosities. 

Another difference emerges when we consider that the lower 
the radio frequency we adopt, the more we are contaminated by the 
lobes which reside further out. To the aim of checking the existence 
and validity of a relation between the jet output at rad io and X-ray 
freque ncies, shorter wavelengths are to be preferred. Panessa et al.l 
( I2OO7I) studied an optically selected sample of ~50 Seyfert galax- 
ies, and compared the 2-10 keV luminosities to the 1.4, 5, and 
15 GHz ones. In particular they found that the best fitting relation, 
taking into account also censored data, is represented by the same 
relation that holds for a sample of LLAGN, with a slope of between 
0.71 and 0.97 at 5 GHz, and 0.8 and 1.02 at 15 GHz, with (without) 
the inclusion of a few radio-loud sources (according to the def- 
inition of Terashima & Wilson 2003) respectively. The difference 
betwe en their Seyfer t samp le and the LLAGN is the mass range, 
which iMerloni et alj ( 120031) showed to be the third parameter in 
the so-called "fundamental plane of black holes" , which is anyhow 
only a scale parameter. It is worth noting that Ide Gasperin et al.l 



8 D. Burlon et al. 



12 






.... 





.... 


J 




p 


■ 




Blazar 


10 


~ 






~ 


8 


- 


/ \ 

1 
1 


- 


Z 6 






1 


- 










4 
2 


/ 


^S 1 


\ 


- / 


r - ^ - 

/ 


\ 





.«-r-rr^. . . 1 . . . . 1 . . . . 1 . ,y,'.'. . . . . 1 rr>r-U-,.i.-« III.. 



-4 -3 -2 -1 

log(Rx) [Rx=L2oghz/Lbat] 

Figure 6. Distributions of tiie Rx parameter. 



( 120 11 1) found only a mild relation in faint LLAGN between 2-10 
keV X-rays and 4.8 GHz and 8.4 GHz. 

The framework of the "fundamental plane of black holes" is 
to establish both theoretically and observationally the coupling be- 
tween the disk (i.e. the accretion power) and the jet. To that aim the 
accretion power linked to the X-ray emission from the corona and 
the radio power is associated to the jet. One of the result of their 
work is to successfully interpret their findings in the framework of 
the fundamental connection between a radiatively inefficient accre- 
tion flow_and_thejeVunde^the assumption of scale invariance of the 
jet JHeinz & Sunvaevll2003 ). From an observational point of view 
they proved that this coupling holds at all scales (i.e. both for AGN 
and Galactic black holes). 

The absence of the radio-X coiTelation in our hard X-ray se- 
lected sample is at odds with the equivalent works at lower radio 
frequencies and lower X-ray energies. In principle this could be 
due to either an incompleteness of our sample (our result arises 
when cross-correlating ~ 20 AGN), or be physical in nature. We 
checked if the relation was weaker or stronger at different radio 
frequencies. Not all sources have counterparts in the AT20G at 
lower frequencies, but as we show in Table |2] the result seems 
to hold at 8 GHz and 5 GHz. We looked at even lower frequen- 
cies, using the fluxe s from the NRAO VLA Sky Survey (NVSS; 
ICond on et al. 1998), Sydney University Molonglo Sky Survey 
(SUMSS; Mauch et al. 2003), and MGPS-2 I Murphy et al. 2007), 
but unfortunately just 8 Seyfert AGN can be securely associated to 
our selection. These numbers are too low to run any meaningful sta- 
tistical analysis. The problem of having too few counterparts to our 
cross-matched sample at lower frequencies, is mirrored at higher 
frequencies. We also checked the 30 GHz and 40 GHz flux densi- 
ties from the Planck satellite Early R elease Compact Source Cata- 
log JPlanck CoUaboration eFZIbOl ih v-l.fl Less than 5 Seyferts 
among our sources have high frequency measurements, making the 
higher frequency analysis even less significant. 

Motivated by the similar result of absence of correlation when 
using the VLBI radio luminosities (Panessa & Giroletti, 2013, MN- 
RAS submitted), we run the following check on our sample. By us - 
ing the definition of compactness adopted by Chhetri et al.J ( l2012h . 
we excluded from our analysis the sources whose visibilities of the 



longest-over-shortest baselines were smaller than 0.86. This effec- 
tively translates in a factor of ~ 10 better angular resolution, down 
to 0.15 arcsec resolutiorjj. Out of the whole sample, 21 AGN sur- 
vive this selection and again the correlation is absent at the same 
statistical level to the original analysis. Therefore we conclude that, 
if the absence of the correlation at high frequencies and high angu- 
lar resolution will be proven statistically solid by future works, it 
could mean that the known radio-X coiTelation is relating a gen- 
uine "core" property (the X-rays) with a more distant one arising 
from the lobes of the jet. This in turn would imply that the X-ray 
output of the AGN, averaged over the duration of the used survey, 
correlate with the energy reservoir represented by the distant radio 
structures. 

Beamed sources: Note that lMerloni et al.l ( 120031) excluded the 
radio-loud (blazar-like) objects, and so do most authors. This is 
because in the blazar case both luminosities are associated to the 
relativistic jet and therefore the electromagnetic output is boosted 
by a factor oc F^. This is not normally the case for Seyfert galax- 
ies, in which the jet -if present- is assumed to be sub-relativistic. 
Testing the universality of the disk-jet coupling at yet other scales 
is beyond the scope of this paper. At the same time, they were mo- 
tivated by the fact that the radio-X connection described by low 
mass systems and high mass systems defined parallel tracks in the 
luminosity plane. Indeed more massive objects like active galaxies 
tend to be more radio-loud with respect to the X-ray output when 
compared to Galactic black holes. The mass segregation and the ac- 
cretion rate are in fact two important parameters of their treatment. 
It is nonetheless interesting to note a similar effect in our sample, 
if we compare the occupation of BAT selected Seyfert-like AGN 
and BAT selected FSRQ. In Fig.[5]we draw with dotted lines differ- 
ent radio-to-X ratios. Seyfert-like AGN tend to populate the area 
where the radio power is between 0.01 and 0.1 the X-ray power, 
while blazars tend to populate the more "radio loud" area where 
the ratio is between 0. 1 and 1. 

Another way to visuali se this effect is to take the definition 
of radio loudness adopted bv lTerashima & Wilsonl (l2003h adopting 
a similar "high frequency" radio loudness Rx = log{Lii/Lx)- 
Fig. |6] shows that the distribution of the Rx parameter is different 
for the two classes of sources, which is a reflection of the two sam- 
ples being displaced along different tracks. The median (and dis- 
persion) of log{Rx) for local Seyferts and blazars are -1.88 (0.69), 
and -0.68 (0.34), respectively. The null hypothesis of the two popu- 
lations of being drawn by the same parent sample has a probability 
P = 1.5 X 10""*. 

Our decision to include also the blazars in this work was to dis- 
cuss the following: if the distribution of the Lorentz factors of the 
jets is peaked (i.e. if there is a "typical" F for BAT FSRQ), then the 
beaming should impact the normalisation of the correlation drawn 
by their luminosities rather then the slope. Conversely, one has to 
consider that probably the dissipation of the jet power takes place in 
different areas along the jet, therefore the beaming will be different 
(higher) for the X-ray luminosities than for the radio luminosities. 
This effect likely affects also the slope of the correlation, but cannot 
be taken into account unless one independently knows the beaming 
factors for radio and X-rays. Even if beyond the goal of this pa- 
per, we would like to mention that if such a control on the Lorentz 
factors could be o btained (e.g. from va riability in X-rays with the 
NuSTAR satellite JHarrison et allbOlOl) , and from VLBI in radio), 



http://irsa.ipac.caltech.edu/Missions/planck.html 



^ This translates into structures more compact than ^^ 9 pc at the median 
Seyfert redshift. and ~ 1.3 kpc at the median blazar redshift 
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then we would be able to de-boost the luminosities of blazars, and 
compare them properly with the Seyfert ones. 



5 SUMMARY AND CONCLUSION 

Starting from the six-year 5wi/f/BAT survey of AGN, we looked 
for counterparts at high frequency radio in the AT20G survey 
JMurphv et alJuOlOl) of the southern sky. By adopting a Bayesan 
cross-matching algorithm we found 37 AGN active in both 
regimes. Both surveys are biased towards bright objects, given that 
the sensitivity of the BAT is of the order of 6 x 10^^^ erg/cm^ /s 
and the AT20G has a cut at 40 mjy. At the same time, the use of 
high energy X-rays is basically unbiased with respect to interven- 
ing gas, while the use of the highest radio frequency available helps 
in the contamination of the radio fluxes by extended structures. This 
is the best available set of observations (excluding VLBI) at these 
frequencies, in order to confront the "core" properties of AGN. We 
summarise our findings in the following: 

• We investigated the ra dio-X connection at -so far- unexplored 
X-ray frequencies (but see Panessa et al.l201ll) : 

• We found that only ~20% (considering that the AT20G covers 
the southern sky) of the hard X-ray AGN have a counterpart at 
20 GHz, and these are equally distributed between local Seyfert 
galaxies and distant blazars; 

• The existence of a relation between the X-ray and radio out- 
put of Seyfert-like AGN in the form Lx oc L%^'^ is essentially 
determined by the redshift; 

• The revised radio loudness parameter shows a bimodality be- 
tween the two classes of sources. A transition between the two 
could be ascribed to the fact that we are using a flux-limited sam- 
ple, beamed luminosities for one class, and to a different mass dis- 
tribution of these sources. 

With the advent of the fi rst focussed hard X-ray mission 
(NuSTAR uHarrison et al.ll2010t) there could be the interesting pos- 
sibility to target the pre-selected sample of Swift/BKT AGN and 
get a lower limit on the Lorentz factor of the high-energy emission. 
From extensive VLBI measurements we woul d be able to resolve 
the ve ry base also of Seyfert radio jets (see e.g. iGiroletti & Panessal 
120091) and even constrain the beaming factor of the radio luminos- 
ity. This will in principle allow the testing of the radio-X connec- 
tion for both Seyfert and blazars at the same time, down to an hard 
X-ray sensitivity 100 times deeper than the current limit. 



ACKNOWLEDGMENTS 

DB is funded through ARC grant DPI 10102034. DB thanks R. Ek- 
ers and F. Panessa for fruitful discussions and useful feedback. This 
research has made use of the NASA/IPAC Extragalactic Database 
(NED) which is operated by the Jet Propulsion Laboratory, Cal- 
ifornia Institute of Technology, under contract with the National 
Aeronautics and Space Administration, and the SIMBAD database, 
operated at CDS, Strasbourg, France. Finally it has benefitted from 
extensive use of TOPCATO 



http://www.star.bristol.ac.uk/~mbt/topcat/ 



APPENDIX A: PARTIAL CORRELATION ANALYSIS 

Given that both Lbat and Lat20 have the same scaling with the 
distance modulus, we tested whethe r the redshift z is driving the 
correlation (see e.g. Padovanilll992l) . as this is a well recognised 
problem when comparing quantities that depend on a third vari- 
able in the same fashion. To this aim we performed a partial cor- 
relation analysis. Renaming the Pearson cortelation coefficients 

^[iBATiAT20] = '^12' "^liBAT^] = "^W' '^[iAT20^1 = '^23, and 

riLoi-rL ^^nn.z] = >'i 2.3 the first order correlation coefficient has 
the form ( iKendall & St uart 1979): 



ri2,3 = 



ri2 — »"i3r23 



To test the significance of a non-zero value for r we compute 



t 



rVJV-2 



which obeys the probability distribution of the Student's t statis- 
tic, with N-2 degrees of freedom, N being the number of elements 
of the samples we are testing. The corresponding probability of a 
chance correlation is computed numerically by integrating the area 
in the two tails of a Student's t-distribution function. The results 
are summarised in Table |2] For each of the three frequencies we 
report in the first and second line the Spearman correlation coef- 
ficient and the corresponding chance correlation probability. In all 
frequencies and subsamples, this statistics alone seems to point to a 
very strong correlation between the observables (as often reported 
in the literature). In the third (labelled with a 2: in parenthesis) and 
fourth line we report the first order Pearson partial correlation co- 
efficient and the coiTesponding probability computed with a Stu- 
dent's t-distribution function (which can be converted into "sig- 
mas" of confidence, should the reader prefer to). 
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